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Abstract 

Where a 2D problem of optimal profile in variable speed flow is 
resolved in a class of convex Bezier curves, using symbolic and numer- 
ical computations. 
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We study a particular problem on searching the optimal 2D shape or pro- 
file of a cartesian object immersed in a variable speed flow and subjected to 
some constraints on the boundary. 

Let be {x, y} a cartesian system, and y = f(x) a function, defined for 
semplicity on the unit interval [0, 1], whose graph is the profile which we want 
to optimize. The function / is subjected to the following constraints, arising 
from engineering requirements: 



cl. /(0) = 
c2. /(1) = 
c3. / is convex in (0, 1) 
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Assume that the freestream flow is a vector field parallel to y-axis with a 
speed distribution according to a function v = v(x) depending on x-variable. 
Also, assume that the flow pressure acting on the object profile is determined 
by Newton's sin-squared law ([2]), that is, being p the fluid constant density 
and a the angle between flow direction and tangent to profile, 

p(x) = pv 2 (x)sin 2 a = pv 2 (x) ] (1) 

1 + f' 2 (x) 

The local force acting on an element profile of x-extension dx is p(x)dx, and 
this is the expression of the local resistance offered to the flow by the body. 
The total force on the profile is therefore 

v 2 (x) 



The purpose of the mathematical analysis is the minimization of this func- 
tional in a suitable class of functions where one hopes to find the solution 
function f(x). In general, usual methods of calculus of variations are applied, 
but one can easily test that technical difficulties arise in simple cases too. For 
example, let be v(x) = ax. The Euler-Lagrange equation (see [3]) associated 
to the minimization problem is than 



n - d dF o 2 d x 2 f 

U " dxdf dx(l + f>2) 2 



from which a four-degree algebraic equations in the variable /' arises, with c 
arbitrary constant: 

* 2 /'+c(i+r) 2 = o (4) 

Such equation is not simple to manage and not so useful. It is more istruc- 
tive, and useful for concrete applications, to consider a class of functions / 
int the domain of Bezier curves. 

Consider the three control points {Pj, % — 0, 1, 2} = {0, 0}, {a, 1},{1, 0}, 
where < a < 1. The associated Bezier curve (see PQ) is 



b(t) = (hit) Mi)) = © ^(l-*) 2 "* = {2at + (l-2a)t 2 ,2t-2t 2 } (5) 

i=0 
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Figure 1: Some Bezier curves in the case a = 0,0.2,0.4,0.6,0.8,1, and relative 
control points. 



Note that the parametric curve b satisfies the conditions cl. and c2., respec- 
tively for t = and t = 1, and c3. (see Fig.(l)). 

Also, x'{t) = 2a + 2(1 — 2a)t, and some usual algebraic computations show 
that x'(t) ^ if 5 < a < 1. Under this assumption, if / = f(x) is a cartesian 
representation of b, we have 



Now the minimization problem is finding the value of the variable a, if exists, 
for which previous functional has its minimum value. 

Consider, for example, the case of a freestream speed distribution v(x) = 
— 5x 3 ; this situation has been proposed by a customer, as velocity field at 
oulet of a ventilating system and going towards a convex body. Along the 
profile of the Bezier curve ([5]) the speed distribution has expression v(x(t)) = 




df db 2 dt db 2 \db{\ 1 l-2t 



(6) 



dx dt dx dt dt a + (1 — 2a)t 

and dx = b[dt, so that the functional ([2]) to minimize becomes 






6 



dt 



(8) 
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We note that previous is an integral of a rational function and it is symbol- 
ically computable, but for our purpose we can only compute it numerically 
for a sufficient set of values of the parameter a, equally spaced in the inter- 
val [0, 1], and verify the existence of a local minimum for F(a). Figure (2) 
shows the evidence of such a minimum for a value 0.6 < ao < 0.8. An exact 
computation, by usual numerical techniques, gives the value ao = 0.682564 
(our computation: use of FindMinimum, by Mathematica v. 7, (see 11])). 



F 




Figure 2: Graph of F(a) in the case v(x) = — 5x 3 . 

The problem of finding the optimal profile in the suitable class of Bezier 
curves is so resolved in this concrete case. 




Figure 3: The optimal profile. The point of maximum is in the region where the 
flow speed values are greater. 
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